Baltic Astronomy, vol. 11, 341-366, 2002. 



KURUCZ MODEL ENERGY DISTRIBUTIONS: 
A COMPARISON WITH REAL STARS. 
II. METAL-DEFICIENT STARS 



(N 
O 

\ V. Straizys 1 , R. Lazauskaite 1 ' 2 and G. Valiauga 1 



1 Institute of Theoretical Physics and Astronomy, Gostauto 12, 
Vilnius 2600, Lithuania; 

2 Department of Theoretical Physics, Vilnius Pedagogical University, 



o 
6 



X 



Studenty, 39, Vilnius 2340, Lithuania. 
Received July 25, 2002. 



> 
CO 

Abstract. Energy distributions of synthetic spectra for Kurucz 
{Nj ■ model atmospheres are compared with observed energy distributions 

of metal-deficient stars of the blue horizontal-branch (BHB), F-G-K 
subdwarf (SD) and G-K giant (MDG) types. The best coincidence 
of the synthetic and observed energy curves is found for BHB stars. 
The largest differences are found in the ultraviolet wavelengths for 
subdwarfs and cool MDGs. The influence of errors of effective tem- 
peratures, gravities and metallicities is estimated. 

Key words: stars: fundamental parameters, spectral energy distri- 
bution - stars: horizontal-branch, subdwarfs, metal-deficient giants 



1. INTRODUCTION 

Synthetic energy distributions of the Kurucz models (Kurucz 
1995) are frequently used for a selection of optimum positions of 
passbands for multicolor photometric systems and for calculation of 
parameters of photometric systems or of transformation equations 
between color indices of different photometric systems. However, the 
use of synthetic spectra for these tasks is justified only in the case 
if the synthetic spectra of model atmospheres coincide sufficiently 
well with energy distributions of real stars having the same physical 
parameters (effective temperatures, gravities and metallicities). 
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In our first paper (Straizys, Liubertas & Valiauga 1997) we have 
compared energy distributions of the Kurucz model atmospheres 
with the mean energy distributions of stars of various spectral types 
of solar metallicity published by Straizys & Sviderskiene (1972) and 
revised by one of the authors. In the present paper we compare the 
Kurucz synthetic spectra of metal-deficient models with spectral en- 
ergy distributions of real stars - blue horizontal-branch (hereafter, 
BHB) stars, F-G subdwarfs (hereafter, SD) and G-K metal-deficient 
giants (hereafter, MDG) published by Sviderskiene (1992). 

2. INTRINSIC ENERGY DISTRIBUTION CURVES OF STARS 

Energy distribution functions of 52 stars - 9 blue horizontal- 
branch stars, 11 extreme F-G-K subdwarfs, 10 mild F-G subdwarfs, 
14 mild metal-deficient G-K giants and subgiants and 6 extreme 
metal-deficient G-K giants - were published by Sviderskiene (1992). 
The wavelength dependence of energy fluxes for each of these stars 
was determined by averaging the data from all sources found in the 
literature until 1990. For the majority of stars the flux data cover 
the range from 300 to 1050 nm. For some stars observed with the 
IUE orbiting observatory, the fluxes are available down to 150-170 
nm. The energy distribution curves are given in the absolute-relative 
form, i.e. they all are normalized to 100 at 550 nm. The mean energy 
distribution of each star is given in the intensities per unit wavelength 
interval. Each value given in the catalog is the intensity integral in 
the interval AA = 5 nm centered on a given wavelength. The stars 
affected by interstellar extinction, are dereddened. The accuracy of 
the observed energy distributions will be estimated in the last section 
of this paper. 

In the catalog of Sviderskiene (1992), Table 4 on pages 44-48 
has a typographic error - the numbers of stars are given in wrong 
positions. The right order of HD numbers in the heading of the 
Table should be: HD 127665, 135722, 216763, 217906, 219615 and 
2857 (instead of 217906, 219615, 2857, 127665, 135722 and 216763). 

3. A COMPARISON OF THE MODEL AND THE OBSERVED 
SPECTRA 

For the present investigation we have selected only the stars 
with extreme metal-deficiency: 6 blue horizontal- branch stars, 11 
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extreme subdwarfs and 6 extreme metal-deficient giants. These stars 
are listed in Table 1. 

Effective temperatures, gravities and metallicities for most of 
these stars are taken from Cenarro et al. (2001) where these pa- 
rameters are reduced to a unique system based on high-dispersion 
spectral analysis. According to that paper, the estimated errors of 
the parameters are ±100 K for T e , ±0.2 dex for log g and ±0.15 
dex for [Fe/H]. For some stars the parameters are taken from other 
sources listed in the Notes to Table 1. 

Sviderskiene (1992) has excluded the interstellar extinction effect 
on energy distribution curves for some giants which were found to 
be slightly reddened. We have made a new search of interstellar 
reddenings for the stars of Table 1. The values of color excesses, 
given in the Table, were used to deredden energy distributions of the 
stars. Intrinsic energy distributions io(A) were calculated from the 
observed / (A) by the equation 

W) = ^± (i) 

where r(A)is the transmittance function of the unit mass of inter- 
stellar dust and x is the amount of unit masses. The values of x for 
different color excesses were calculated by the equation 

X =^Y ( 2 ) 

(x=l) 

Here E B _ V ' for stars of various spectral types are calculated from 
the standard interstellar extinction law r(A) given in Table 3 of the 
Straizys (1992) monograph using the following equations: 

_ fl(\)R m (\)T*(\)d\ 

2 - 51 ° g fl(\)R m (\)d\ ■ ^ 

and 

-^7711— 7712 Ani ^-777,2 • (4) 

Here / (A) are the intrinsic energy distribution functions of stars, 
-R m (A) are the response functions of the B and V passbands, mi 
and 1712 are B and V magnitudes. Mean energy distribution func- 
tions for stars of various spectral classes were taken from Straizys & 
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Sviderskiene (1972) with the slightly revised ultraviolet. The values 

(x=l) 

of E B _ V due to the band-width effect slightly depend on the intrin- 
sic energy distribution of the star: for O-B-A stars (including BHB 
stars) it is close to 1.00 and for the G5-K5 giants (including MDG 
stars) it is ~0.85. 

For the stars, which were dereddened by Sviderskiene, only dif- 
ferential reddening due to different color excess value was taken into 
account. 

Figures 1-23 compare the energy distribution curves of Kurucz 
and Sviderskiene for the metal-deficient stars of Table 1. The tem- 
peratures of the models are rounded to the nearest number multiple 
of 10. For each star two Kurucz energy curves are plotted with the 
temperatures ±200 K from the right value. The space between these 
two curves is shown in grey. The synthetic curves for the specific 
T e , log g and [Fe/H] are obtained by linear interpolation between the 
grid curves given by Kurucz. The following conclusions can be made 
for different types of stars. 

(1) For the BHB stars the coincidence of energy distributions is 
relatively good. However some stars show considerable differences in 
some spectral ranges. For example, the model of HD 60778 (Figure 2) 
is stronger by 10-20% in the ultraviolet and violet and fainter by 
~5% in the red and infrared. This difference can be corrected by 
taking a lower value of T e of the star. The model energy curve of 
HD 161817 (Figure 6) is also higher by ~10% in the ultraviolet and 
violet, but no difference is seen in long wavelengths. 

(2) In case of the extreme subdwarfs the coincidence is relatively 
good in the whole spectral range (from the ultraviolet to the infrared) 
only for HD 25329 (Figure 8) and BD+26 2606 (Figure 16). For the 
remaining nine subdwarfs the model curves, corresponding to the T e 
values given in Table 1, are too low by 10-15% in the blue, violet 
and ultraviolet ranges. These differences cannot be explained by a 
T e error, since in the red and infrared the model and the observed 
curves are close to each other. Only for HD 19445 (Figure 7) an 
increase in the temperature of the star by 200 K would force the 
observed and the model curves to be in sufficient agreement. 

(3) For the extreme metal-deficient giants the coincidence of the 
observed and model curves in the red and infrared is also good. For 
HD 2665 (Figure 18) and HD 165195 (Figure 22) the model curve is 
good also in the violet and ultraviolet. However, for the other four 
stars there are systematic differences in the short wavelength part of 
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Table 1. A list of the stars whose spectral energy distribution data are 
compared with the Kurucz models. 



HD, BD 




l°g 9 


[Fp / Hi 
[re/nj 


Eb— V 


3HB stars 










2857 


7563 


2.0" 


-1.60 


0.02 


60778" 


8690 


3.30 


-0.50 


0.01 


74721 b 


8640 


3.55 


1 AO 

-1.48 


0.01 


86986 b 


7850 


3.10 


-1.81 


0.02 


109995 


8034 


2.98 


-1.55 


0.00 


161817 


7639 


2.96 


-0.95 


0.01 


F-G-K subdwarfs 








19445 


5918 


4.35 


-2.05 


0.00 


25329 


4787 


4.58 


-1.72 


0.00 


64090 c 


5446 


A A P* 

4.45 


-1.7b 


0.00 


84937 


6228 


A ( \ 1 

4.01 


-2.17 


0.00 


94028 


5941 


A 01 

4.21 


1 /in 

-1.49 


0.00 


103095 


5025 


4.56 


-1.36 


0.00 


140283 


5687 


3.55 


-2.53 


0.01 


188510 


5490 


4.69 


-1.59 


0.00 


219617 


5878 


4.04 


-1.39 


0.00 


+26 2606 d 


6146 


4.23 


-2.29 


0.01 


+17 4708 


6005 


4.01 


-1.74 


0.01 


Gl-K metal-deficient giants 








2665 


5013 


2.35 


-1.96 


0.06 


6755 


5102 


2.40 


-1.41 


0.04 


88609 


4513 


1.26 


-2.64 


0.02 


122563 


4566 


1.12 


-2.63 


0.01 


165195 


4471 


1.11 


-2.15 


0.15 


221170 c 


4465 


1.04 


-2.10 


0.03 



NOTES: 

° HD 60778: parameters are from Danford & Lea (1981). 
6 HD 74721 and HD 86986: parameters are from Adelman & Philip (1994). 
c HD 64090 and HD 221170: parameters are from Soubiran et al. (1998). 
d BD +26 2606: parameters are from Gratton et al. (1996). 
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the spectrum: there is a general weakening of model spectra in the 
violet and ultraviolet for HD 6755, HD 88609, HD 122563 and HD 
221170. The upper limit of the shadowed area in short wavelengths 
almost coincides with the real energy curve for all these four stars. 
This might indicate that T e accepted for the stars is by ~200 K 
too low. However, acceptance of higher temperature will lead to a 
disagreement of energy curves in the red and infrared. 

For the evaluation of the effect of gravity uncertainties, we have 
chosen three stars, representing BHB, SD and MDG stars. The value 
of logg was varied by ±0.3 dex from the true value, keeping the T e 
and [Fe/H] values constant. Model energy curves for these stars are 
compared with the observed energy curves in Figures 24-26. It is 
evident that the gravity effects are of minor importance and they 
cannot help to explain the differences seen in Figs. 1-23. The same 
is true in the case of [Fe/H] variations by ±0.3 dex (Figures 27-29). 

4. A COMPARISON OF SYNTHETIC COLOR INDICES 

For a quantitative evaluation of the differences between the 
model and the observed energy distributions, we calculated and in- 
tercompared synthetic color indices in the Vilnius photometric sys- 
tem and the UBV system of the metal-deficient stars from Table 1 
and the stars of luminosities V and III of solar chemical composition 
(hereafter, normal stars). The method of synthetic photometry used 
is described by Straizys (1996). Energy distributions for normal stars 
were taken from Straizys & Sviderskiene (1972) (with the corrected 
ultraviolet, as described in Paper I) and for metal-deficient stars - 
from Sviderskiene (1992). The response functions of the Vilnius and 
UBV systems were taken from Straizys (1992). 

The results are shown in Figs. 30-37, where the differences 
A of corresponding color indices ("model" minus "observed") are 
plotted against the effective temperatures. The temperatures for 
metal-deficient stars are taken from Table 1 and for the normal stars 
- from Paper I. 

Synthetic color indices reflect the differences seen in energy dis- 
tributions. For stars of solar chemical composition (dots and crosses) 
the scatter of stars gradually increases with decreasing temperature. 
However, for stars with T e >4500 K the deviations A do not exceed 
±0.2 mag even in case of color indices containing the ultraviolet 
magnitudes. The deviations of subdwarfs are also within these lim- 
its. The largest deviations (plus 0.3-0.4 mag) are observed for the 
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three coolest metal-deficient giants with the temperatures near 4500 
K. It is probable, that in the ultraviolet these stars are really brighter 
than their Kurucz models. 

5. DISCUSSION AND CONCLUSIONS 

The comparison of energy distributions of real stars and their 
models, described in the previous sections, is a difficult problem for 
interpretation since there are several causes affecting the accuracy 
of the results. They are: (1) accidental and systematic errors of the 
observed energy distributions, (2) systematic errors in the account 
of opacity sources in model stellar atmospheres, and (3) errors of 
the accepted effective temperatures of the stars. Let us discuss these 
error sources separately. 

When we plot energy distributions of a star, determined by dif- 
ferent authors, and normalize all of them at one wavelength, say at 
550 nm, then we find a considerable scatter of individual determi- 
nations at both ends of the wavelength interval covered. If we limit 
ourselves with the traditional ground-based data, covering the wave- 
lengths between 320 nm and 1 /j,m, then the usual scatter of the 
energy flux values for G-K-M type stars is: ±5-10% at 350 nm and 
±5% at 1 (Straizys & Sviderskiene 1972, Sviderskiene 1988). The 
flux values obtained as averages from the results of several authors 
should be more accurate. 

A good test of the accuracy of the average spectral energy curves 
is a comparison of color indices obtained from these curves by syn- 
thetic photometry and color indices of the same stars observed pho- 
toelectrically with sufficiently high accuracy (±0.01 mag). Such a 
comparison in the passbands of the Vilnius photometric system has 
been done both for normal stars of various luminosities (Sviderskiene 
1988) and for metal-deficient stars (Sviderskiene 1992). In both cases 
the stars exhibit differences between the ultraviolet synthetic and ob- 
served color indices within ±0.10 mag. Consequently, a considerable 
part of the differences seen in Figures 1-37 can be explained by the 
errors in observed energy distributions. 

It is strange enough that the differences of the observed and 
the model energy distributions give an evidence that at least part of 
metal-deficient stars (subdwarfs and giants) exhibit overestimation 
of opacity in stellar atmospheres. It is more usual that synthetic 
spectra of model atmospheres demonstrate a "missing opacity" . On 
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the other hand, some model stars exhibit a perfect coincidence of 
their spectral energy distributions with observations. If the opacity 
were a problem, then all stars of a given type would behave in the 
same way. 

The errors of the accepted effective temperatures is the third 
possibility. Although Cenarro et al. (2001) estimate their tempera- 
tures to be accurate within ±100 K, it is not excluded that in some 
cases these errors may be twice larger. The vertical bars in Figures 
35-37 show how the temperature errors of ±200 K affect the plotted 
differences of color indices. It seems that these error bars are capa- 
ble to explain the scatter of most stars in these figures, except of the 
stars cooler than ~4500 K. 

Thus we conclude that the present investigation gives evidence 
that spectral energy distributions of the Kurucz metal-deficient mod- 
els of the effective temperatures >4500 K do not contradict the 
observed energy distributions of real stars belonging to the blue 
horizontal branch, extreme subdwarfs and extreme metal-deficient 
giants, if we take into account the errors of the observational data 
and of the effective temperatures. Cooler models probably exhibit 
some systematic differences in the ultraviolet and violet wavelengths. 
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Fig. 1. Spectral energy distribution of HD 2857 (BHB star). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 2. Spectral energy distribution of HD 60778 (BHB star). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 3. Spectral energy distribution of HD 74721 (BHB star). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 4. Spectral energy distribution of HD 86986 (BHB star). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 5. Spectral energy distribution of HD 109995 (BHB star). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 6. Spectral energy distribution of HD 161817 (BHB star). The 
effective temperature of the model is changed by ±200 K. 



352 



V. Straizys, R. Lazauskaite, G. Valiauga 





Fig. 8. Spectral energy distribution of HD 25329 (subdwarf). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 9. Spectral energy distribution of HD 64090 (subdwarf). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 10. Spectral energy distribution of HD 84937 (subdwarf). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 12. Spectral energy distribution of HD 103095 (subdwarf). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 14. Spectral energy distribution of HD 188510 (subdwarf). The 
effective temperature of the model is changed by ±200 K. 



356 



V. Straizys, R. Lazauskaite, G. Valiauga 




300 400 500 600 700 800 

X(nm) 

Fig. 15. Spectral energy distribution of HD 219617 (subdwarf). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 16. Spectral energy distribution of BD+26 2606 (subdwarf). 
The effective temperature of the model is changed by ±200 K. 
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Fig. 17. Spectral energy distribution of BD+17 4708 (subdwarf). 
The effective temperature of the model is changed by ±200 K. 
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Fig. 18. Spectral energy distribution of HD 2665 (MD giant). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 20. Spectral energy distribution of HD 88609 (MD giant). The 
effective temperature of the model is changed by ±200 K. 
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Fig. 22. Spectral energy distribution of HD 165195 (MD giant). 
The effective temperature of the model is changed by ±200 K. 
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300 400 500 600 700 800 

l(nm) 

Fig. 23. Spectral energy distribution of HD 221170 (MD giant). 
The effective temperature of the model is changed by ±200 K. 




Fig. 24. Spectral energy distribution of HD 109995 (BHB star). 
The gravity of the model is changed by ±0.3 dex. 
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Fig. 26. Spectral energy distribution of HD 165195 (MD giant). 
The gravity of the model is changed by ±0.3 dex. 
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Fig. 27. Spectral energy distribution of HD 109995 (BHB star). 
The metallicity of the model is changed by ±0.3 dex. 
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Fig. 28. Spectral energy distribution of BD+26 2606 (subdwarf). 
The metallicity of the model is changed by ±0.3 dex. 
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Fig. 29. Spectral energy distribution of HD 165195 (MD giant). 
The metallicity of the model is changed by ±0.3 dex. 
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Fig. 30. A comparison of the U—V color indices of stars and their 
models (model minus observed). Mean wavelengths of the U and V pass- 
bands are 345 and 544 nm. The vertical bars indicate the errors originating 
from the temperature changes by ±200 K. 
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Fig. 31. The same as in Fig. 30 but for P-V. The mean wavelength 
of the P passband is 375 nm. 
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Fig. 32. The same as in Fig. 30 but for X—V. The mean wavelength 
of the X passband is 405 nm. 
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Fig. 33. The same as in Fig. 30 but for Y—V. The mean wavelength 
of the Y passband is 466 nm. 
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Fig. 34. The same as in Fig. 30 but for Z— V. The mean wavelength 
of the Z passband is 516 nm. 
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Fig. 35. The same as in Fig. 30 but for V-S. The mean wavelength 
of the S passband is 656 nm. 
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Fig. 36. The same as in Fig. 30 but for U-V of the UBV system. 
Mean wavelengths of the U and V passbands are 364 and 550 nm. 
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Fig. 37. The same as in Fig. 30 but for B-V oi the UBV system. 
The mean wavelength of the B passband is 442 nm. 



